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£N) ' Abstract 



The Painleve first equation can be represented as the equation of 
/"*\ . isomonodromic deformation of a Schrodinger equation with a cubic 

potential. We introduce a new algorithm for computing the direct 
monodromy problem for this Schrodinger equation. The algorithm is 
based on the geometric theory of Schrodinger equation due to Nevan- 
linna 



< 

u 

43* 1 Introduction 

-i— > 

a . 

Painleve equations form the core of what may be called "modern special 
function theory" . Indeed, since the 1970s' pioneering works of Ablowitz 
and Segur, and McCoy, Tracy and Wu the Painleve functions have 

C^) • been playing the same role in nonlinear mathematical physics that the 

classical special functions, such as Airy functions, Bessel functions, 
etc., are playing in linear physics. 

ly-s , For example, it has been recently conjectured by Dubrovin |Dub08| 

(and proven in a particular case by Claeys and Grava |CG09j) that 

I — [ , Painleve equations play a big role also in the theory of nonlinear waves 

f~^ ' and dispersive equations. 

^^ , In this context, the first Painleve equation (P-I) 

y" = 6y 2 -z, zeC , (1) 

- 1 — i i 
S^ ■ is of special importance. 

5h \ Indeed recently [DGK09J Dubrovin, Grava and Klein discovered 

that a special solution of P-I -called integrate tritronquee- provides 
the universal correction to the semiclassical limit of solutions to the 
focusing nonlinear Schrodinger equation. 

The key mathematical fact about the Painleve equations is their 
Lax-integrability. This allows to apply to the study of the Painleve 
functions a powerful Isomonodromy - Riemann-Hilbert method. Using 
this method, a great number of analytical, and especially asymptotic 
results have been obtained during the last two decades- for the general 
theory see the monumental book [FIKN06 , for what concerns P-I see 
our papers (MaslOa [M aslObj . 
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It is clear that we need efficient and reliable algorithms for comput- 
ing solutions of the Painleve equations and that these algorithms should 
take into account the integrability of Painleve equations. To this aim 
S. Olver has recently build an algorithm to solve the Riemann-Hilbert 
problem -or inverse monodromy problem- associated to Painleve-II 
[OlvlOj . 

In the present paper we introduce a new simple algorithm for solv- 
ing the direct monodromy problem associated to P-I; given the Cauchy 
data y,y',z of a given solution of P-I we can compute its monodromy 
data -or equivalently, the Riemann-Hilbert problem associated to it. 

Given our algorithm, it would be possible to compute the inverse 
monodromy problem for Painleve-I simply using an appropriate New- 
ton's method. Due to the simplicity of our algorithm, this procedure 
would be probably more efficient than a numerical solution of the corre- 
sponding Riemann-Hilbert problem (which, to the best of our knowl- 
edge, has not yet been given). However, this is a matter of future 
research. 

Below we briefly introduce the direct monodromy problem associ- 
ated to P-I. 

1.1 The Perturbed Cubic Oscillator 

Consider the following Schrodinger equation with a cubic potential 
(plus a fuchsian singularity) 

^^ = Q(X;y,y',zMX), (2) 

Q(X;y,y',z) = 4A 3 - 2Az + 2zy - Ay 3 + y 



a , v' 



X-y 4(A-y) 2 



We call such equation the perturbed (cubic) oscillator. 
We define the Stokes Sector Sk as 



s k = U 



27rfc 

arg A — 

5 



< J |> , k e z 5 ■ 



Here, and for the rest of the paper, Z5 is the group of the integers 
modulo five. We will often choose as representatives of Z5 the numbers 
-2,-1,0,1,2. 

For any Stokes sector, there is a unique (up to a multiplicative 
constant) solution of the perturbed oscillator that decays exponentially 
inside Sk- We call such solution the k-th subdominant solution and let 
i/) k (X;y,y',z) denote it 0. 

The asymptotic behaviour of ipk is known explicitely in a bigger 



1 Equation ([2|) has a fuchsian singularity at A = y. Hence, as we explain in Section 2, 
any solution is two-valued having a square root singularity at A = y. We do not discuss 
this fact here: the reader can find the rigorous definition of subdominant solution in the 
Appendix. 



sector of the complex plane, namely Sk-i U Sk U Sk+i- 

lim M^f) _>i >Ve >0. (4) 

| argA ^|< ¥ _ £ A-iexp{-|Al + fA^} 

Here the branch of A 2 is chosen such that tpk 1S exponentially small in 

s k . 

Since ipk— 1 grows exponentially in Sk, then ifik—i and ipk are linearly 
independent. Then {ip^-i, tpk} is a basis of solutions, whose asymptotic 
behaviours is known in Sk-i U Sk- 

Fixed fe* £ Z5, we know the asymptotic behaviour of {ipk*-i,'4 , k»} 
only in Sk*-i U Sk* ■ If we want to know the asymptotic behaviours of 
this basis in all the complex plane, it is sufficient to know the linear 
transformation from basis {ipk-i,ipk} to basis {ipk,ipk+i} for any k £ 
Z 5 . 

From the asymptotic behaviours, it follows that these changes of 
basis are triangular matrices: for any k, tpk—i — V-'fe+i + < J fcV'fc for some 
complex number o~k , called Stokes multiplier. The quintuplet of Stokes 
multipliers o~k , k £ Z5 is called the monodromy data of the perturbed 
oscillator. 

We can now define the direct monodromy problem. 

Direct Monodromy Problem. Fixed y, y', z compute the Stokes 
multipliers ak(y, J/', z) of the perturbed oscillator equation (0). 

Our Algorithm gives a numerical solution of this problem. 

1.2 P-I: Isomonodromy Approach 

What is the relation of the direct monodromy problem with P-I? 
Painleve-I is the equation of isomonodromic deformation of 
the perturbed oscillator. Indeed the following Theorem holds. 

Theorem 1. Let the parameters y = y(z),y' — v } z ' of the potential 
Q(\; y, y' , z) be functions of z; then y(z) solves P-I if and only if the 
Stokes multipliers of the perturbed oscillator do not depend on z. 

Proof. See for example [MaslOa . □ 

Fix a solution y(z) of P-I. If z is a pole of y(z) then equation (|5} 
is not well-defined. However, recently the author [MaslOa (sec also 
|CC94| ) showed that this difficulty can be overcome. 

Let a £ C be a pole of y(z) then y(z) has the following Laurent 
expansion around a 

, x 1 a(z — a] 2 (z — a) 3 ,. , 4 _,, .*. 

2/0) = + y — ' + - — —!- + b{z - a) 4 + 0((z - a) 5 ) . 



(z-a) 2 ' 10 6 

Inserting the Laurent expansion into equation @ and taking the 
the limit z — > a, equation ^ becomes the following Schrodinger equa- 
tion with a cubic potential 

d f( A ^ = V(X; 2a, 286WA) , V(X; a, b) = 4A 3 - aX - b . (5) 
dX z 



Moreover the Stokes multipliers of the equation (J5J) coincide with the 
Stokes multipliers of equation ([2]). 

We call equation ([5]) the cubic oscillator. For the rest of the paper 
we consider it as a particular case of the perturbed oscillator ||2J) . 

1.3 Stokes Multipliers and Asymptotic Values 

The algorithm is based on the formula © below, that we discovered 
in [MaslOcj . 

Consider the following Schwarzian equation 

{f(X),X} = -2Q(\;y,y',z). (6) 

Here {/(A), A} = pA/ — § ( pA) ) is the Schwarzian derivative. 
For every solution of the Schwarzian equation ([5]) the following limit 

exists 

w k (f)= lim /(A)eCUoo, 

A — >oo ,A£Sfc 

provided the limit is taken along a curve non-tangential to the bound- 
ary of S k . 

In Section 2, we will prove that the following formula holds for any 
solution of the Schwarzian equation ([5]) 



vk{y,y',z) = i(wi + k(f),w-2+k(f);w-i+k(f),w2+k(f)) ■ (7) 

(a-c){b-d) 
(a-d)(b-c) 



Here (a, b; c, d) = (_ < 7\/ h _ I is the cross ratio of four points on the 



sphere. 

The paper is organized as follows. In Section 2 we introduce the 
Nevanlinna's theory of the cubic oscillator and the Schwarzian differ- 
ential equation ^ . Then we prove formula ([7]) formula for computing 
the Stokes multipliers from any solution of the Schwarzian differential 
equation. Section 3 is devoted to the description of the Algorithm. In 
Section 4 we test our algorithm against the WKB prediction and the 
Deformed TBA equations. For convenience of the reader, we explain 
the basic theory of cubic oscillators (Stokes sectors, Stokes multipliers, 
subdominant solutions, etc ...) in the Appendix. 
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2 Schwarzian Differential Equation 

As we mentioned in the Introduction, our Algorithm is based on for- 
mula (J7J| that allows one to compute Stokes multipliers from any solu- 
tion of the Schwarzian differential equation ([5]) . 

This formula, proven in Theorem [5] below, has its roots in the ge- 
ometric theory of the Schrodinger equation, which was developed by 
Nevanlinna in the 1930s' |Nev32| . The author learned such a beautiful 
theory from the remarkable paper of Eremenko and Gabrielov |EG09| . 
In this section we follow quite closely |EG09| as well as author's recent 
paper |MaslOc| . 

REMARK. Equation 0) has a fuchsian singularity at the pole A = y 
of the potential Q(X;y,y' , z). However this is an apparent singularity 
JMaslOaf : the monodromy around the singularity of any solution of 
(0) is —1. As a consequence, the ratio of two solutions of (0) is a 
meromorphic function. 

The main geometric object of Nevanlinna's theory is the Schwarzian 
derivative of a (non constant) meromorphic function /(A) 

if(\)\\- f " {X) 3 r /,/(A) V (8) 

{/(A) ' A} -/W~H7W • (8) 

The Schwarzian derivative is strictly related to the Schrodinger 
equation ([S]). Indeed, the following Lemma is true. 

Lemma 1. The (non constant) meromorphic function f : C — > C 
solves the Schwarzian differential equation 

{f(X),X} = -2Q(X;y,y',z). (9) 

iff /(A) = ^j-j4 where 0(A) and x(A) are two linearly independent so- 
lutions of the Schrodinger equation fU). Hence, the first derivative of 
any (non constant) solution of |[P|) vanishes only at the pole A = y of 
the potential. 

We define the Asymptotic Stokes Sector Sk as 

2irk 



Sk = A : 



argA- 

5 



< Ji> 



Every solution of the Schwarzian equation ^ has limit for A — > oo, 
A G Sk- More precisely we have the following 

Lemma 2 (Nevanlinna). (i) Let /(A) = -^Wy be a solution of 0) 
then for all k G Z5 the following limit exists 

w k (f)= lim /(A) eCU 00, (11) 

A — >-oo A£S/b 

provided the limit is taken along a curve non-tangential to the 
boundary of Sk ■ 

(a) wk+i{f) ^ w k (f) , Vfc e Z 5 . 



(Hi) Let g(X) = f$$ = ggggg , (^ J) € 01(2, C). T/>en 

«*(*) = ^#^ • (12) 

CWkif) + « 

('toJ J/ £/ie function f is evaluated along a ray contained in Sk, the 
convergence to Wk(f) is super-exponential. 

Proof, (i-iii) Let ipk be the solution of equation @ subdominant in Sk 
and i/>fc+i be the one subdominant in Sk+i- Since they form a ba- 

sis of solutions, then /(A) = ^gl±g^W > for some f « ^ 
GJ(2,C). Hence io fc (/) = f if S ^ 0, w k {f) = oo if 5 = 0. Sim- 
ilarly tifc+i(/) = f • Since (° ^ G Gl(2,C) then Wfc (/) ^ 

wfc+i(/) 

(iv) From estimates ((3]) we know that inside S^, 

ipk(X;y,y',z) 



ipk+i(\;y,y',z) 



-Ke( |A5- 2 A2 



where the branch of A 2 is chosen such that the exponential is 
decaying. 

□ 

Definition 1. Let /(A) be a solution of the Schwarzian equation (0) 
and Wk{f) be defined as in ill}) . We call Wk(f) the k-th asymptotic 
value of f . 

We noticed in a previous paper [Mas 10c| that the Stokes multipliers 
of the Schrodinger equation are rational functions of the asymptotic 
values Wk(f)- This relation is the basis of our Algorithm. 

Theorem 2. J.MaslOc'j Denote au the k-th Stokes multiplier of the 
Schrodinger equation H|) (for its precise definition, see equation H8\) 
in the Appendix). Let f be any solution of the Schwarzian equation 
@). Then 

(T k =i(w 1+k (f),w^2+k(f);w-i+k(f) 1 w 2 +k(f)) ,VfceZ 5 , (13) 

where (a, b; c, d) = (a-d)(b-c) * s ^he cross ratio of four point on the 
sphere. 



Proof. Due to equation (fl2|) all the asymptotic values of two different 
solutions of §§§ are related by the same fractional linear transformation. 
As it is well-known, the cross ratios of four points of the sphere is 
invariant if all the points are transformed by the same fractional linear 
transformation. Hence the right-hand side of (|13|) does not depend on 
the choice of the solution of the Schwarzian equation. 

Let ipk+i be the solution of ^ subdominant in Sk+i and 4>k+2 
be the one subdominant in Sk+2 (see the Appendix for the precise 



G 



definition). By choosing /(A) = ^ fc+1 )J , one verifies easily that the 
identity (JT2J) is satisfied. □ 

REMARK. The same construction presented here holds for (inhar- 
monic oscillators with polynomial potentials of any degree. For any 
degree, there are formulas similar to H3\) for expressing Stokes multi- 
pliers in terms of cross ratios of asymptotic values. The general for- 
mula will be given in a subsequent publication. 

2.1 Singularities 

Since the Schwarzian differential equation is linearized (see Lemma [1} 
by the Schrodinger equation, any solution is a meromorphic function 
and has an infinite number of poles (for a proof of this fact see |Nev32) 
and [Elf34j). The poles, however, are localized near the boundaries of 
the Stokes sectors Sk, k € Z5. Indeed, using the estimates (jl]) one can 
prove the following 

Lemma 3. Let /(A) be any solution of the Schwarzian equation |Pp. 
Fix e > and define S k — {A : |arg A - 2|£| < § - e} , k e Z 5 . Then 
/(A) has a finite number of poles inside Sk- Hence, there are a finite 
number of rays inside Sk on which /(A) has a singularity. 



3 The Algorithm 

In the previous section we have proved the following remarkable facts 

• Along any ray contained in the Stokes Sector Sk, any solution 
/ to the Schwarzian differential equation (J5|) / converges super- 
exponentially to the asymptotic value Wk{f)- See Lemma |5](iv). 

• The Stokes multipliers of the Schrodinger equation © are cross- 
ratios of the asymptotic values Wk{f)- See equation (|13[) . 

• Inside any closed subsector of Sk, f has a finite number of poles. 
See Lemma [3] 

Hence the Simple Algorithm for Computing Stokes Multipliers goes 
as follows: 

1. Set k=-2. 

2. Fix arbitrary Cauchy data of /: /(A*), /'(A*), /"(A*), with the 
conditions A* ^ y, /'(A*) ^ 0. 

3. Choose an angle a inside Sk, such that the singular point X = y 
does not belong to the corresponding ray, i.e. a 7^ argy. Define 
t : R+ U -> C, t(x) = f(e ia x + A*). The function t satisfies the 
following Cauchy problem 

f {t(x),x} = e 2ia Q(e ia x + \*;y,y',z), 

1 t(0) = /(A*), i'(0) = e M /'(A*), i"(0) = e 2ia f"(\*) . 



4. Integrate equation ([T4")) either directly [j or by linearization (see 
Remark below), and compute Wk(f) with the desired accuracy 
and precision. 

5. If k < 2, k + +, return to point 3. 

6. Compute ai using formula (fl~3|) for all I £ Z5. 



REMARK. As was shown in Lemma{^ any solution f of the Schwarzian 
equation is the ratio of two solutions of the Schrodinger equation. Hence, 
one can solve the nonlinear Cauchy problem |i^| ) D U solving two linear 
Cauchy problems. 

Whether the linearization is more efficient than the direct integra- 
tion of \1$ will not be investigated in the present paper. 



4 A Test 

We have implemented our algorithm using MATHEMATICA's ODE 
solver NDSOLVE integrating equation (fl"4"|) with steps of length 0.1. 
We decided the integrator to stop at step n if 

|t(0.1n)-t(0.1(n-l))| < 10~ 13 and 



i(0.1n)-t(0.1(n- 1)) 



i(0.1n) 



< 10 -13 



To test our algorithm we computed the Stokes multiplier ao(b) of 
the equation 

^ = (4A 3 - bMX) . (15) 

According to the WKB analysis (see |Sib75| . |Masl0a| ) the Stokes 
multiplier 00(b) has the following asymptotics 



Ti-r(i/3) 5 
be 



_j e ?7rr(ii/6)"" if b > 



1 • • (-")" { v/fr(i/3) , -\ 

cos { 2^r(n/6) { - b) ~ 6 )> iib<0 - 

(16) 
Our computations (see Figure 1 and 2 below) shows clearly that 
the WKB approximation is very efficient also for small value of the 
parameter b. 

We also tested our results against the numerical solution (due to A. 
Moro and the author) of the Deformed Thermodynamic Bethe Ansatz 
equations (Deformed TBA), which has been recently introduced by 
the author [MaslOcJ, developing the seminal work of Dorey and Tateo 
|DT99| . The Deformed TBA equations are a set of nonlinear integral 



integrating equation H14|l directly, one can hit a singularity x* of y. To continue the 
ition past the pole, starting from x* — e one ca 
satisfies the same Schwarzian differential equation. 



solution past the pole, starting from x* — e one can integrate the function y — -, which 
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Figure 1: Thick dotted line: the rescaled Stokes multiplier i e 23r ( 11 /e) ao(b) eval- 



uated with our algorithm; thin continuous line: cos 
prediction for the rescaled Stokes multiplier. 
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Figure 2: Thick dotted line: the rescaled Stokes multiplier e 22/3r C 11 /« cr (&) eval- 
uated with our algorithm; thin continuous line: 1, the WKB prediction for the rescaled 
Stokes multiplier. 



equations which describe the exact correction to the WKB asymptotics. 
The numerical solution of the Deformed TBA equations enable to a- 
priori set the absolute error in the evaluation of the Stokes multiplier 
(To (6) rescaled with respect to the exponential factor, shown in (flTj)) . 
Hence, in the range of —20 < b < 20 we could verify that we had 
computed the rescaled ao{b) with an absolute error less than 10 -8 . 

5 Appendix 

The reader expert in anharmonic oscillators theory will skip this Ap- 
pendix; for her, it will be enough to know that we denote (Tk{y,y',z) 
the k — th Stokes multipliers of equation @. Here we review briefly 
the standard way, i.e. by means of Stokes multipliers, of introducing 
the monodromy problem for equation ([2|) . All the statements of this 
section are proved in Appendix A of author's paper [MaslOaj and in 
Sibuya's book |Sib75| . 

Lemma 4. Fix k G Z5 = {—2, . . . , 2} and define a cut in the C plane 
connecting A = y with infinity such that its points eventually do not 
belong to S k —i U S k U Sk+i- Choose the branch of X? by requiring 

lim RcA^ = +00 , 

A— >oo 
argA=2|i 

while choose arbitrarily one of the branch 0/A 3 . Then there exists a 
unique solution if)f.(X;a,b) of equation |3J) such that 

lim MX -ff Z \ ^l,Ve>0. (17) 

Definition. We denote ^ k the k-th subdominant solution or the solu- 
tion subdominant in the k-th sector. 

From the asymptotics (jTTJ) , it follows that ip k and ipk+i are linearly 
independent. If one fixes the same branch of A * in the asymptotics 
(JTTJ) of ipk—ij V'fej V'fe+l then the following equations hold true 

ipk-i(My,y',z) ipk+i{^;y,y',z) , , , 

il>k{X;y,y',z) tp k (X;y,y',z) 

(18) 
-icr k+3 (y,y',z) = 1 + <r k (y, y' , z)cr k+1 (y, y' , z) , Vfc £ Z 5 . 

Definition. The entire functions a k (y,y' , z) are called Stokes multi- 
pliers. The quintuplet of Stokes multipliers a k (y,y' , z),k £ Z5 is called 
the monodromy data of equation (J5j). 
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